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Abstract 

The problem of P vs. NP is very serious, and 
solutions to the problem can help save lives. This 
article is an attempt at solving the problem using 
a computer algorithm. It is presented in a fashion 
that will hopefully allow for easy understanding 
for many people and scientists from many diverse 
fields. 

In technical terms, a novel method for solving 
k-SAT is explained. This method is primarily 
based on linear algebra and finite fields. Ev- 
idence is given that this method may require 
rougly 0{ii?) time and space for deterministic 
models. More specifically the algorithm runs in 
time 0{P ■ V{n + VY) with mistaking satisfiable 
Boolean expressions as unsatisfiable with an ap- 
proximate probablity l/<d{V{n + VY)^ , where 
n is the number of clauses and V is the num- 
ber of variables. It's concluded that significant 
evidence exists that P=NP. 

There is a forum devoted to this paper at 
http://482527.ForumRomanum.com. All are in- 
vited to correspond here and help with the 
analysis of the algorithm. Source code for 
the associated algorithm can be found at 
https://sourceforge.net/p/la3sat. 



1 Introduction 

There are many problems proposed to computer sci- 
entists that have been thought to be too difficult for 
computers to solve quickly. In fact, perhaps the most 



fundamental question in computer science is to find if 
certain types of problems, collectively known as the 
class NP, can be solved quickly by a computer. If so, 
a world of opportunities would open up, and many 
new problems that were supposed to be almost im- 
possible to solve could be solved quickly. This paper 
attempts to provide a proof that they can be solved 
quickly, and also shows a way to do it. This will 
hopefully invite researchers from many diverse fields 
to contribute to the research and work of solving NP 
hard problems. 

The level of interest in this question is so great that 
in 2000, the Clay Mathematics Institute listed the 7 
Millennium Prize Problems, and offered $1,000,000 
for someone who could prove the relationship between 
P and NP f7] , which would answer the question. 

1.1 Turing Machines 

To understand the classes P and NP, first consider 
the basic notion of a Turing machine, which is a sci- 
entific definition of a computer. It has one or more 
tapes that hold information. At any time, the Turing 
machine uses the tape or tapes to decide what to do 
next. The machine has a prescribed set of instruc- 
tions to it to help it decide. 

To develop the ideas behind Turing machines more, 
a distinction is made between what types of deci- 
sions can be made. A deterministic Turing machine 
(DTM) has a predetermined decision for every type of 
situation it encounters; thus the term deterministic. 
A nondeterministic Turing machine (NTM), on the 
other hand, can have more than one action it can do 
for any given situation. In other words, it's actions 
aren't determined ahead of time. Therefore, it's said 
that it's nondeterminstic. 

Figure [T] makes the difference more clear. NTMs 
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A decision tree is shown below. As an algorithm (computer program) progresses, it must pick a 
choice from three possibilities. Nondeterministic turing machines(NTMs) have been thought to have 
an advantage in this case, because, at each step, they can pick from any choice. Deterministic turing 
machines (DTMs) don't really have the ability to pick, so they're thought to be at a disadvantage. 



START 



O Choice available to NTM 

O Choice available to NTM and DTM 

Figure 1: A Scenario Involving Choices 



have been thought to be able to solve more problems 
than DTMs because of the availability to make more 
choices. Here the DTM has only one choice avail- 
able at each step, while the DTM has three choices 
available. So this evidence points to DTMs as be- 
ing more limited in potential than NTMs. In fact, 
there is enough difference between these two com- 
puters that different classes of computation (compu- 
tational power) were proposed for each. 



1.2 P and NP 

Right around 1970, Leonid Levin, in the USSR, and 
Stephen Cook, of the US, independently arrived at 
the concept of NP-completenessp]. 

The idea of NP-completeness comes from two 
classes of algorithms. Once again, these two classes 
come from the ideas of DTMs and NTMs. 

The classes are defined partially by how much can 
be accomplished within a limited amount of time. 
This time limit, which will be explained briefly here, 
is defined by the mathematics of asymptotic analy- 
sis, which is described in fUJ . Basically, the time 
limit is given as a function of the problem size. The 
function can take on many forms, and two common 
forms are shown above in Figure[2j Here a polynomial 
function (t = 0{n°')) is contrasted with an exponen- 
tial function {t = 0(a")). Note that the polynomial 
functions may take more time for smaller problems, 
but the larger problems, which computer scientists 
are mainly concerned with, have smaller times for 
polynomial functions. 

The two classes mentioned above are defined for 
polynomial functions. One class, P, is defined as all 
problems that can be solved by a DTM in polynomial 
time (a polynomial function of the problem size). The 
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The difference between functions 



other class, NP, is defined as the class of all prob- 
lems that can be solved by an NTM in polynomial 
time. A^P-complete problems are then problems in 
NP that are at least as hard to solve as the hardest 
problems in NP. 

The fundamental question that this paper attempts 
to answer is if A^P-complete problems can be solved 
in polynomial time by a DTM; in other words, is 
P = NPl 

For more information on P versus NP, one can 
refer to Sipser ^6j. 

1.3 SAT and k-SAT 

SAT essentially asks if, given a Boolean formula, can 
the formula be satisfied. In other words, it asks if the 
variables in the formula be given a truth assignment 
that makes the entire formula true. It's required, to 
thoroughly answer the question, that a certificate is 
given if the answer is true. This is usually in the 
form of a particular solution to the question, such as 
a satisfying assignment for all variables. 

k-SAT is a particular variation of SAT, in which the 
formula is organized into clauses. All variables inside 
a clause are connected via disjunction. All clauses 
are connected via conjunction. The k in k-SAT then 
refers to the number of variables in each clause. More 
information on SAT and k-SAT can be found in [T^] . 

Conventional SAT and k-SAT solvers function by 
learning clauses or making random guesses at solu- 



tions ([14] and the recent survey [15]). 



The current best upper bounds for 3-SAT is 
0(1.32113") time, and can be found in ISAAC 2010 
by Iwama et al[TB]. There's an arXived paper by 
Hertli, Moser and Scheder which gives an 0(1.321") 
time algorithm [17]. These are all randomized algo- 
rithms. There is an arXived, deterministic 3-SAT 
algorithm by Kutzkov and Scheder that runs in time 
0(1.439")[B|. 



There is a paper that proposed a polynomial run- 
time for a SAT variant. V.F. Romanov proposed a 3- 
SAT algorithm that uses "discordant structures" , or 
set-like operations on a lattice, to represent informa- 
tion and determine solutions in polynomial time|19). 
However, Liao presents an argument that P is not 
equal to NP, using a 3-SAT variant, in [9]. 



Perhaps Baker, Gill, and Solovay's theory of rela- 
tivization helps explain why there are not more algo- 
rithms that attempt to solve NP-complete problems 
in polynomial time [20]. Their paper "Relativizations 
of the P != NP Question" states that consulting an 
oracle can lead to situations in which P != NP. 
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Figure 3: Representations of xi 



2 Data Structure(s) 

One of the most fundamental components of the al- 
gorithm will be refered to as the clause polynomial. 
Figure [3] shows the basic organization of clause poly- 
nomials. Essentially, the complete information of one 
clause is represented by a polynomial of one vari- 
able. That is, for eaxh particular truth assignment 
of all variables, the polynomial "records" whether the 
clause satisfies this assignment. This is then attached 
to the polynomial's variable, which orders the truth 
assignments. 

Note that the particular truth assingnments are 
shown in red. In this fashion, the same tree that 
organizes the variables can be repeatedly used, and 
therefore provides some standard organization to the 
information. Again, this tree organizes all possible 
variable assignements. It does so by starting at the 
root, and then proceeds through all Boolean variables 
in order, and assigns one node to each possible assign- 
ment of true or false for each variable. The leaves (at 
the bottom) thus represent a complete truth assing- 
ment for all variables. 

The clause polynomial for a single variable is equiv- 
alent to a digit in a binary number. In order to under- 
stand this. Figure |4] shows binary numbers composed 
of one to three bits. The tree in Figure |3] uses a three 
bit example, since there are three variables; thus one 



1 Bit 


2 Bits 


3 Bits 


0(0) 


00(0) 


oOo(o) 


1(1) 


01(1) 


oOi(i) 




10(2) 


0l0(2) 




11(3) 


0ll(3) 






l00(4) 






lOl(5) 






ll0(6) 






lll(7) 



Figure 4: Binary Representations 
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xo In System Of xq, xi, X2 






x^\l+x^'){l+x'") = 






x^ + x^ + x^ + x"^ = 






0x0 ^ l^.! ^ 1^3 ^ 0x4 _^ 1^5 ^ 0x6 ^ 1^7 








xi In System Of Xq, xi, X2 






(l+x2°)x2'(l+x2') = 






x2 + X^ + X^ + x"^ = 






0x0 ^ 0x1 ^ 1^2 ^ 1^3 ^ 0x4 _,_ 0a;5 ^ ^ I2.7 




Figure 5: Constructing Polynomials For Variables 



bit for each variable. Note then, that the middle vari- 
able, xi, out of xo, xi, and X2, corresponds with the 
middle bit. This can be seen in the figure by looking 
at the middle bit for each binary representation of 
three bits (which is displayed larger than the other 
two bits). The sequence is the same as the binary 
tree above. 

2.1 Constructing Clause 
Polynomials 

A basic clause is of the form 

Xao VXai V---VXa, (2.1) 

The algorithm seeks to construct clause polynomials 
for clauses in this form. This is a two step process. 
The basic theory behind the process is fairly easily ex- 
plained here, although the technicalities and a proof 
are saved for Appendix |X] on page [26l 

To begin to understand how clause polynomials are 
constructed, a few patterns are observed. Figure |6l 
on the next page, helps to demonstrate these. First, 
the obvious pattern of binary digits appears for the 
variables on the lefthand side. All variables have a 
pattern of 2' zeros and then 2* ones for variable Xi, 
which then repeats. These variables, in conjunction, 
correspond with binary numbers. In other words, for 
the least significant bit, the pattern is zero, then one, 
and then it repeats. For the next least significant bit 
(or variable), the pattern has two zeros and then two 
ones, which then repeats. This pattern continues for 
all variables. 



Observing this pattern, an easy way to construct 
the clause polynomials for a single variable becomes 
clear. In fact, it can be summarized as a simple for- 
mula: 



\*:=0 / \k=m+l / 

(2.2) 

Here Y[ {f{x)) denotes the (indefinite) product of a 
function of x. More information can be found about 
indefinite products in [26 . This is taken over a sys- 
tem of n -|- 1 variables. 

Figure [S] shows how the clause polynomials are ac- 
tually constructed using the formula. Note that these 
are all for a single variable. For the xq example, 
m = since the variable is xq. In other words, if 
the variable was X32, then the algorithm would plug 
in m = 32 into the formula. Note, then, that n is one 
less than the total number of variables, so in both 
cases it is two, since there are three variables. Re- 
turning to the Xo example, the lefthand product has 
nothing inside it, the middle follows from knowing 
n, and the right side product is determined from m 
and n and the rules of products. For the xi example, 
there is a product on the left side to work with, since 
m is now one, and k is set to zero for this. Then it 
follows that (1 -f x2'°) = (1-1- x^) since fc 0. The 
right side follows similarly, this time with fc — 2. 

Returning attention to Figure [6] on the previous 
page, a second pattern can be observed. This is 
appearant in the blue portion of the figure. Here, 
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Figure 6: Repeating Patterns of Clauses 



For Clauses of the Form V a:„, V • • • V x, 



For each Xa, 

Calculate g{xa,) = + x^")^ x^"' 

Let Result = 
For h = Q to z 

If {h = ai) FOR SOME 

Result = Result + g{xh) 

Else 

Result = Result- (1 + x^") 
Return Result 



Figure 7: Clause Polynomial Construction Algorithm 



Towards P^NP via k-SAT 



©Matt GrofF 2011 



7 





(xo V Xi) From xq, Xi, X2 




g{xo) = -- 


= 1x1 g{xi) = (1 + x^°^ x^' = Ix^ ^ 


r Ix^ 


f{xo V xi) = 






f{xo V xi) = 


Oa;° + Ix^ + la;2 + Ix^ + Ox* + Ix^ + la;^ - 


^Ix' 



{xq V Xi V X2) From xq, Xi, X2 

ff(xo) = = 1x1 ^ 5(xi)= (l+x2'')x2^ =1x2 + 1x3 

g{x2) = (1 + 3^^") (1 + X^^) = 12^" + 12;^ + l2^^ + 12;^ 

/(xo V xi V X2) = (xi) + (x^ + x^) + (x^ + x^ + x^ + x^) 

/(xo V xi V X2) = Ox° + 1x1 + Ix^ + Ix^ + Ix^ + Ix^ + Ix^ + Ix'^ 



Figure 8: Clause Polynomial Examples 



the second pattern is shown in boxes for two- variable 
clauses. There is the original one- variable pattern of 
ones and zeros, which is then followed by a series of 
ones. Note that the series of ones is exactly the same 
size as the original single-variable pattern. This en- 
tire pattern then repeats. 

The observed pattern comes from the combination 
of two variables. The long series of ones comes from 
the variable that is "larger" than the other. The 
short series of repeating ones and zeros comes from 
the "smaller" of the two variables. Together, they 
form a series that repeats. Another way of looking at 
it is that both variables form a repeating series, and 
combining these two repeating series creates another 
repeating series. 

In fact, the pattern of repetition continues even as 
more variables are added to the clause. The purple 
portion of the figure shows the repitition(s) involved 
with three variables. Again, combining two variables 
produces a short, repeating series. When a third vari- 
able (that is strictly "larger" than the other two) is 
added, a longer, yet repeating, series emerges. 

Again, the techinicalities of all of this are presented 
(and solved) in Appendix]^ on page[26l 

There is a very simple algorithm to calculate clause 
polynomials for any amount of variables, as long as 
they are in the form given in Equation 12.11 Figure 
[7] shows the algorithm. The basic idea is to first cal- 
culate the the (possibly long) sequence of ones that 



repeats for each variable. This is then shifted into 
the proper position. It is identical to the clause cal- 
culations for one variable, with the exception that 
the series is not made to repeat. The idea is that us- 
ing "smaller" variables, the algorithm will construct 
short repeating sequences, and add in the appropri- 
ate larger variable when the repeating sequence gets 
large enough. So essentially, it constructs one large 
repeating sequence by making the smaller variable se- 
quences repeat, and adding in larger variables when 
the sequence gets large enough. 



Figure |5] shows examples of creating clause polyno- 
mials for given problems. Note that the complete set 
of variables for the original problem must be known 
ahead of time, just as for clauses of individual vari- 
ables. As can be seen, the results from this figure 
correspond with the first eight values from the corre- 
sponding clauses in Figure [6] 



There is really only one operation remaining to be- 
ing able to construct essentially any clause polyno- 
mial of the form in Equation 12.11 negation. As it 
turns out, negation is not much more difficult than 
constructing non- negated clause polynomials. 
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Figure 9: Patterns of Clauses With Negation 



For Clauses of the Form V V • 




5ra = (nro'(i+^''")) 




If {h = at) for some al 




Result = Result- + gix^) 




Else If {h — a^) for some 




Result = Result + g{xh) 




Else 




Result = Result- (1 + x^'^) 





Figure 10: Negated Clause Polynomial Construction Algorithm 
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(xo V Xi) From Xn, Xi,X2 






g{xo) = = 1x1 g(^) + 2,2"^ ^ ^ 12,1 






/(xo V Xi) = (((lxl)x2') + (1x0 + 1x1)) • (l + X2') 






/(xo V Xi) = + 1x1 + 0x2 ^ 1^3 _^ 1^5 0a;6 _^_ 








[xq V Xi V X2) From xq, Xi, X2 






g{x^) = = 1x1 5(3^) = (1 + a;2") = Ix" + 1x1 
g(x2) = (1 + a;^") (1 + a;^') a;^' = la;'* + la;^ + la;^ + la;"^ 






/(xo V xi V X2) = (^((lxi)x2') + (lx° + 1x1)) • (1 + x2') + (Ix^ + Ix^ + Ix^ + Ix^) 






/(xo V xi V X2) = Ix" + 1x1 + Ox^ + Ix^ + Ix'' + Ix^ + Ix^ + Ix"^ 




Figure 11: Clause Polynomial Examples With Negation 



Figure m on the previous page, displays clause pat- 
terns with negated variables. As can be noted by 
the negated variable Xq near the left, it is just a re- 
versed version of ones and zeros, compared to the 
non-negated xq displayed just to the left of it. In 
other words, just as negation reverses true and false 
values in Boolean algebra, negation reverses zeros and 
ones in clause polynomials. In fact, this is how all 
single-variable clauses are negated; simply swap the 
zeros and ones. There is actually a very simple and 
effective way to do this - the negated clauses are the 
same as the regular clauses except for the fact that 
they are not multiplied by a power of x^ , as the 
regular clauses are. This is the equivalent as remov- 
ing this term from Equation 12. 2[ hence the resulting 
equation: 

(m— 1 \ / " \ 

11(1 + ^'') n + (2-3) 
k=a ) \fe=m+l / 

Again, this is for a system of n-|- 1 variables, and n 
can be adjusted accordingly for the size of the prob- 
lem. 

Next, attention can be turned to the new pat- 
terns observed for multiple-variable clauses. Figure 
[9] shows the negated values of xq next to the com- 
bined clause Xo V xi. Note that the pattern of ones 
coming from xi has not changed in any way; however. 



the repeating pattern of zeros and ones coming from 
Xo has been switched, compared to the non-negated 
equation xq Vxi beside it on the right. So really, it is 
becoming evident that negating variables only swaps 
the pattern of ones that the variable creates. 

This is more evident if the clause Xo V xT is ob- 
served beside it. Here the pattern of ones coming 
from the variable x\ has changed, since x\ has been 
negated. So again, there is evidence that the pattern 
of ones coming from a variable is simply "exchanged" 
or placed elsewhere. It can be seen, though, that this 
does not change the pattern for X2, as evidenced in 
the equation xo V xT V X2 beside it. 

Finally, it's seen that swapping the pattern of ones 
may entail moving a smaller group of ones and zeros. 
This is evidenced in the final clause in the figure. 
Here, the variable X2 has been negated, which causes 
a shift in the largest pattern of ones, which is due to 

X2. 

Two examples have been worked out in slight detail 
in Figure 1111 Here the pattern is mostly the same, 
with the exception that negating variables causes 
some values to be swapped. This is the equivalent 
of moving the pattern of ones. The adjustments to 
the algorithm are shown in Figure [TOl on the previous 
page. 
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Figure 12: Two Clause Satisfaction 



3 Two Clause Problems 

Now that sufficient background has been presented, 
the major ideas behind the algorithm can be intro- 
duced. Two clause problems are some of the simplest 
cases, yet they allow the fundamental ideas to be in- 
troduced and studied. 

Figure [T^] shows the four basic possible combina- 
tions between two clauses. If a clause polynomial is 
considered with any amount of terms (powers of x), 
there are only two possible values (known in math- 
ematics as coefficients) that can be associated with 
each term; zero or one. If two clauses are considered, 
there are two possible values for the coefficient of the 
first polynomial, and two possible values for the co- 
efficient of the second polynomial, for a total of four 
different combinations. 

The figure shows two clause polynomials and the 
trees associated with them. In the middle, the four 
possible combinations can be seen. Remembering 
that the coefficients represent satisfaction for a par- 
ticular truth assignment, it's possible to note the 
total satisfaction for two clauses considered simul- 
taneously. When both clauses have a zero coeffi- 
cient, that indicates that neither clause will satisfy 
the original question of satisfaction for that partic- 
ular truth assignment. Looking through the cluase 
trees, it's seen that this assignment corresponds to 
{xq = false, xi — false). So this truth assignment 
won't satisfy any of the clauses. On the other hand, 
the two truth assignments in the middle of the fig- 
ure each satisfy one clause. The truth assignment 



(xo = true,a;i — true) is seen to satisfy both clauses, 
and so it is a solution to the original problem. Again, 
the reason it satisfies both clauses is that it has a one 
for both coefficients, thus signifying that both clauses 
are satisfied. 

One important thing to note here is that between 
two clauses, there are really only three possible sat- 
isfaction results. Either zero, one, or both of the 
clauses are satisfied (for any particular truth assign- 
ment). One general idea that, although perhaps triv- 
ial, will be important is that there are really only 
three types of satisfaction between two clauses. This 
basic concept will be extended to see that there are 
really only n+1 possible types of satisfaction between 
n clauses. 

The fundamental idea here is that the problems 
can be simplified so that large problems can really be 
dealt with by simply considering the types of satis- 
faction, which are fairly simple. For two clause prob- 
lems, it will really only be necessary to deal with the 
three types of satisfaction, and interactions between 
the different types of satisfaction. 

It will be shown how operations of multiplication 
and addition can be used to work with the three dif- 
ferent types of satisfaction, and to solve questions 
about them. This will be essentially simpler and in 
some ways necessary to overcome the complications 
of working with all of the possibilities between two or 
more clauses. 
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Ox" + Ix^ + 0^2 + Ix^ 



f{xo) = Ox" + 1x1 -^^3 

/(1)-(1+X2'')(1+X2^) 
= (1 + X)(1 +x2) 

= lx° + Ix^ + Ix^ + Ix^ 

h{xo)^{a-f{xo}) + {fil)-f{xo)) 
= (Ox° + ax^ + Ox^ + ax^) + 

(lx° + Ox^ + Ix^ + Ox^) 
= (lx° + ax^ + Ix^ + ax^) 



Figure 13: Example Pre-Multiplication Calculation 



3.1 Manipulating Clause 
Polynomials 

In order to simplify clause calculations into the the 
equivalence classes of satisfaction, it is necessary to 
modify the original clause polynomials. 

The first modification used is a simple one. The 
algorithm will need to change all of the coefficients 
that are equal to one into coefficients that are equal 
to a. This is easy; it is just multiplication by a. Here's 
an example: 

/(x) ^ Ox" + Ix^ + Ox^ + Ix^ 
a ■ /(x) = Ox° + ax^ + Ox^ + ax^ 

The next modification is a bit tougher; it is ex- 
changing the one and zero coefficients. This can be 
done by subtracting the original function from a func- 
tion of ones. This inverts the bits, since the ones are 
subtracted from ones to become zeros, and the the 
zeros are subtracted from ones to become ones. How- 
ever, the function of all ones is needed. This function 
is, for a system of V variables: 

/(i)= na+x^'") (3.1) 

A;=0 

Note that: 

/(I) = Ix" + Ix^ + Ix^ -I- • • ■ Ix^'^-i 

This completes the requirements for multiplication. 
The algorithm uses this knowledge to change the one 



coefficients into a coefficients and the zero coefficients 
into one coefficients. Then multiplication can be per- 
formed, as will be seen. 



FigurelTSlexhibits the steps taken before a function 
is ready for multiplication, /(xq) is given through a 
tree, before multiplication. The ones must be trans- 
formed into as, and the zeros must be transformed 
into ones. Note the final result, h{xo). It is the fin- 
ished calculation, ready for multiplication. 



To do this, the function of ones for a two- variable 
system must be calculated. Then the original func- 
tion can be multiplied by a to transform the one co- 
efficients, and it is also subtracted from the function 
of ones (seperately) to transform the zero coefficients. 
These are added together for the final result. 



Note that this procedure works for any clause in 
any system with any amount of variables. It simply 
prepares the clause polynomials for special process- 
ing, or interactions, with other clause polynomials. 
This may not seem like much, but it greatly simpli- 
fies the system. 
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Xo In System Of xq, xi, X2 

[x^"f (l + {x^"f) (l + {x^'f) = 
x"^ + x^ + x^^ + x" = 

_^ 1^2 ^ ^ 1^6 ^ ^ I^-IO ^ 0x12 ^ 1^14 

(xo V V X2) From xq, Xi, X2 

g(xo) = (x2")2 = Ix^ 5(xi) = (1 + (x2")2) (x^^)^ = Ix^ + Ix^ 

g(x2) = (1 + (x2")2) (1 + (x2')2) (x2')2 = Ix^ + Ix^O + Ix^^ + Ix^^ 

/(xo V xi V X2) = (x^) + (x"* + x^) + (x^ + x^° + x^^ + a;i4) 

/(xo V xi V X2) = Ox° + Ix^ + 1x4 + Ix^ + Ix** + Ix^O + 1x^2 + 1x^4 

Figure 14: Example Pre- Addition Calculation 



Addition requires modifications too. However, 
these modifications are of a different variety than 
those of multiplication. Essentially, the power of x 
needs to be doubled. 

Doubling the power of x can actually be fairly sim- 
ple. Figure [14] displays a retake of the creation of 
clause polynomials. The first part of the example is 
really just Figure [5] on pageO modified for doubling 
the power of x. Note that each power of x, once it is 
essentially calculated, is doubled. The result is that 
all of the powers of x are multiplied by two in the 
finished clause polynomial. 

The second part of the figure is also a redo, this 
time of a multiple variable clause taken from Figure 
[8] on page [71 Again, the main idea is that the powers 
of the variable x is doubled. 



v-i 

filmodified) = n (1 + i^^"?) (3-2) 

k=0 



The actual operation of addition is also performed 
on one special value, which comes in part from the 
ideas of multiplication. A constant is added to each 
coefficient of the polynomials. This constant is the 
same value for all coefficients, so the function of ones 
once again becomes useful. Unfortunately, in its orig- 
inally derived form, the powers of x are not the same 
as the powers of x used in addition. So once again, 
the same principles are used to double the powers of 
X for the function of ones. It is really as simple to 
do this as altering every power of x in the original 

equation (Equation [B1 on page [26]). The new formula These operations will be very useful in the material 
for the modified function of ones is as follows: ahead. 
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Figure 15: Clause Multiplication 
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Figure 16: Operation Equivalence 



3.2 Multiplication 

It could be said that the algorithm is focused around 
(arithmetic) multiplication of clause polynomials. 

Figure [TS] details an example. It begins with two 
clause polynomials, f{xo) and f{xi). The associated 
trees are shown that go with the clause polynomials. 
In the middle, the clause polynomials are broken into 
pieces and multiplied piecewise. Every part of each 
polynomial is multiplied with every part of the other 
polynomial. This corresponds with plain old arith- 
metic multiplication of two polynomials. However, 
an interesting event happens here. The result that 
lies along the diagonal (shown in yellow) corresponds 
with the result f{xo A xi). That is, the diagonal rep- 
resents the corresponding clause polynomial for the 
result of conjunction. The idea here is that multipli- 
cation of clause polynomials can be used to evaluate 
the original Boolean equation. 

Figure [16] helps give a partial explanation of why 
this occurs. As can be seen from the truth tables, the 



logical operation of conjunction (A) and the arith- 
metic operation of multiplication(-) are equivalent 
on bit values. So it's not totally unpredictable that 
multiplication of a clause polynomial contains results 
similar to conjunction. The diagonal in Figure \W\ 
contains like terms multiplied by like terms. Only 
the coefficients (or numbers attached to the terms) 
may differ. The result (along the diagonal) is known 
in mathematics as the Hadamard product, and the 
algorithm is concentrated on seperating this diagonal 
term from the off-diagonal terms. 



The reason why the algorithm is so closely asso- 
ciated with the Hadamard product, or the diagonal 
terms, is that it is essentially the solution to the origi- 
nal problem. Since it tells which terms are satisfying, 
the algorithm can simply count the number of sat- 
isfying terms along the diagonal. If there are any 
satisfying terms, then the original problem can be 
satisfied. Otherwise, it can't be satisfied. 



So at this point in the ideaology, the original 
Boolean problem has been transformed from a ques- 
tion about Boolean arithmetic, to a question concern- 
ing how to get information about the diagonal (or 
Hadamard product). 
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1-x a- X 
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Figure 17: Multiplication and Satisfaction 



3.3 Multiplication and 
Satisfaction 

Figure [17] shows another viewpoint of arithmetic mul- 
tiplication. This time the four different cases based 
on the coefficients are shown. In other words, be- 
tween any two coefficients being multiplied, the input 
coefficients must be one of four cases. Also, on the 
left, the satisfaction between the two original coeffi- 
cients is shown. Then, on the right, the correspond- 
ing modified coefficients and their results are shown. 
Note that there are still really three cases of satisfac- 
tion, resulting in Ix^, ax"^ , and a^x'^. 

It's seen that modifying the coefficients for multi- 
plication has allowed for the three cases to appear in 
the results of multiplication. This is one of the keys 
to getting things to work correctly. The algorithm 
is really interested in the case when both clauses are 
satisfied. 

Unfortunately, the results for multiplication mix 
the diagonal and off-diagonal cases together. In or- 
der to isolate the diagonal, there will be some tricky 
interplay between multiplication and addition ahead. 

One thing that can be noted is that multiplica- 
tion with the original clause polynomials isolates the 
case where everything is maximally satisfied (both 
clauses are satisfied). This can be seen in Figure 
[THl Note that only the case on the far right (where 
both clauses have nonzero coefficients) returns any- 
thing other than zero. Now this will occur for both 




diagonal and off-diagonal values, but it's very close 
to a solution. The algorithm wants the diagonal por- 
tion of this result, and seeks to somehow eliminate 
the off-diagonal portion of it. 



The next subsection will show how the algo- 
rithm can begin to seperate diagonal terms from 
off-diagonal terms. The algorithm's efforts will be 
concentrated on seperating diagonal values from off- 
diagonal values, since the result will lead to a solu- 
tion. 
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co(lx , aQX , aox , ao x ) 


EE [cox'',coaox^,coao^x^) 


ci{lx^, aix^, oix^, ai^x^) 


— [cix ,ciaix ,ciai X ) 


(Oa;2,la;2,la;2,2x2) 




+ 


= {dx'^,{d+l)x'^,{d + 2)x^) 


{dx^,dx'^,dx'^,dx'^) 






(0, 0, 0) + off-diagonal 



Figure 19: Combining Operations For Elimination 



3.4 Eliminating The Diagonal 

Figure [20l on the following page, presents the cases 
for addition alongside the better-explored operation 
of multiplication. The main thing to note is that ad- 
dition can return one of three results, just like the 
equivalence classes of satisfaction. In fact, they are 
once again related in the same way that multiplica- 
tion is related to satisfaction. 

Figure [H] shows what is proposed for the algo- 
rithm. It will take two multiplication operations 
(three are shown, but this is more than needed), and 
combine them with one addition, to eliminate ev- 
erything except for the ofF-dagonal values, which are 
shown as gray triangles. 

Everything is really summarized as arithmetic in 
Figure [191 Here, the original results are shown on 
the left, being modified so that they can be com- 
bined together. It can be noted that the middle two 
cases have always had equal results, so they have been 
combined. The right side shows only three cases in 
paranthesis, which are the satisfaction equivalences. 
The algorithm seeks to eliminate these, thus the sums 
at the bottom are zero, except for the off-diagonal. 

Here the algorithm comes up with three equations 
that must be satisfied: 

co + ci-d = (3.3) 
coao -1- cifli - (d -Fl) = (3.4) 
coao^-f cifli^ - (d + 2) =0 (3.5) 

Again, these come from the three cases on the right 
side of Figure 1191 summing the columns. 

To summarize all of this again, what is essentially 
happening is that two multiplications, together with 
an addition, cancel out all coefficients along the di- 
agonal. However, multiplication creates off-diagonal 



values, and these values will remain. So in effect, the 
algorithm isolates off-diagonal values. It does so by 
using the values calculated from addition to cancel 
out the diagonal from multiplication. 

So now the algorithm can concentrate fully on elim- 
inating the terms along the diagonal to isolate the off- 
diagonal terms. The following subsection will explore 
how to get these terms to cancel. 

First, it's time to introduce modular arithmetic. 
Subsection IC.llI in the Appendix notes some sources 
that go over modular arithmetic. To simplify the 
ideas, essentially all calculations are performed as 
usual, except that an additional operation is per- 
formed afterwards. After the calculations are done, 
to get the result modulo a prime p, the algorithm 
does the equivalent of taking the remainder after di- 
viding the result by p. Thusly, all calculations are 
represented by an integer greater than or equal to 
zero and less than p. So all calculations have a very 
limited range of values that can result. 

This introduces a notion of equivalence, where two 
values are equivalent if they are the same modulo 
p. This allows the algorithm to find solutions more 
easily, since the normal restriction that calculations 
must be equal is relaxed so that calculations must 
only be equivalent. 

Equivalence is really introduced so that the equa- 
tions that must be satisfied can be solved. The 
relaxed restrictions allow for solutions that can be 
found easily. 
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Figure 21: Eliminating The Diagonal 
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Figure 22: Multiplication Results Modulo 7 



To find solutions that eliminate the diagonal, a fi- 
nite field is introduced, allowing operations to be con- 
ducted modulo a prime p. Figure [21] displays calcu- 
lations inside a field modulo 7. Here the algorithm 
selects ao = 2, although it could select really any 
value other than one or zero. It then calculates the 
results of multiplication (times a constant cq), which 
for the various equivalence classes are Co • 1, Co • ag, 
and Co • ao^. 

The algorithm is really interested in the second- 
order difference. This is illustrated in Figure [531 
That is, it takes the difference between results, and 
then takes the difference of these differences. 

It's fairly straightfoward to get the initial results; 
using arithmetic multiplication works fine. Then the 
results have to be modulated. As an example, in 
Figure [231 ai is set as three. Then, to calculate ai^, 
take 3^ = 9. Then 9/2 = 1 with remainder 2. So the 
result is the remainder, which is 2. 

That gets the initial results, in the field. Then to 
get the first order differences, take successive results 
and subtract the first from the successive. This is il- 
lustrated very clearly in the figure, which does a bet- 
ter job of explaining. The take the difference again, 
which is called the second-order difference. 

The reason that it's so important to get these 
second-order differences is that they help match up 
multiplication results. Note that both figures are ac- 
tually two seperate multiplications; Figure lo uses ao 
and Co while Figure [23] uses ai and ci. They both use 
the same prime, so they are two seperate multiplica- 
tion results that can be matched up. The goal is to 













1 = 1 mod 7 1 3 _ ^ 

ai = 3 mod 7 , 

1 2 — 3 = 
ai^ = 2 mod 7 ' 


= 2 , 

} 6 - 2 = 4 
= 6 ' 




Figure 23: DifTerences 
Results 


of Multiplication 



find the second-order differences that add up to zero 
mod p. Note that the example with ai has a second- 
order difference of four. 3 + 4 = 7 = mod 7, so the 
idea is to find a second-order difference of three in 
the other equation. It is done with co = 3. So now 
two equations have been found that match up. 

To check this result, the equations can be com- 
bined: 

co-l-Hci-l = 3- l + l- l =4 = d + 0e 
Co • ao + ci • ai = 3 • 2 + 1 • 3 = 2 = d + le 
Co • oo^ ci • ai^ = 3 • 4 1 ■ 2 =0 = d + 2e 

The main thing to note is that the sequence of re- 
sults in the equations (4, 2, 0) can be recreated by an 
addition operation (on clause polynomials). That is, 
the first equation = 4-f-0(-2), the second = 4+l(-2), 
and the third = 4-f 2(-2). 
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Figure 24: Two Clause Results 



Figure [23] shows the results of combining equations 
(addition and multiplication). Here the the results 
of clause operations (multiplication and addition) are 
combined, and are shown according to whether or not 
they lie along the diagonal. As mentioned previously, 
the addition operation cancels out the multiplication 
along the diagonal. However, the results off the diag- 
onal are not cancelled. 

Thusly, the results off the diagonal become multi- 
pliers for unknown quantities. This is because the 
polynomials that underline the equations are also 
part of this mix, and the resulting polynomial quan- 
tities are unknown, even though the multipliers are 
known. Thus these quantities are represented as bo, 
6i, and 62. 

This allows for a new result. Remember that the 
original quantities are combined together in the previ- 
ous equations that were used. Thus, this new result 
represents a sum of the new quantities, along with 
their respective multipliers. Therefore, the algorithm 
arrives at a new result which is a single equation in 
three unknowns. 

Again, time for some ideas. Back in Section 13.31 
on page[T4l Figure ITSl is discussed. Especially impor- 
tant is the fact that the maximally satisfied clauses, 
both along the diagonal and off the diagonal, can be 
isolated. In this section it has been shown that equa- 
tions for off-diagonal values can be obtained. Think- 
ing about off-diagonal values, it is possible to cre- 



ate other equations using multiplication and addi- 
tion - and these equations can help determine the 
off-diagonal values by combining them with the first 
equation and then using linear algebra to determine 
the values. 



So essentially, the algorithm will use multiplication 
and addition of modified clause polynomials to cre- 
ate equations in three unknowns. These unknowns 
are the off-diagonal values. Then, the equations cre- 
ated can be combined via linear algebra to detetmine 
the off-diagonal values. As mentioned previously, the 
algorithm can already isolate the maximally satis- 
fied clause values together. Unfortunately, the off- 
diagonal and diagonal values are combined together. 
However, since the off-diagonal portion can be deter- 
mined via linear algebra, the diagonal portion can be 
isolated from the off-diagonal portion using simple al- 
gebra, and the results from linear algebra (which give 
the off-diagonal portion). 



The diagonal value of maximally satisfied clauses is 
thus determined, and this corresponds to the number 
of satisfied solutions. Now if this number is nonzero, 
the whole problem can be satisfied. Otherwise, the 
original problem can't be satisfied. 
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4 A Two Clause 
Walkthrough 

In order to help understand the algorithm, this walk- 
through will start at the beginning of execution and 
proceed through as many steps as possible. The ini- 
tial problem will be given as the equation: 

(xo) A (x^Vxi) (4.1) 

One of the first things to do is to determine the 
modulus. Using a prime greater than (2n)^ should 
work well, and it's obvious that there are 2 clauses 
and 2 variables. 17 should suffice, since 17 > (2 • 2)^. 

Next, a value should be given to x, which will be 
used to calculate the clauses. There doesn't seem 
to be any particularly good way to pick, other than 
using a number greater than one. Three will be used 
for this example. 

Next, the values for the clauses can be calculated. 
Recalling Equation 12.21 on page [5j 

/(-™)= fnV+-^^))-^" f n 

\fc=0 / \k=m+l ) 

(4.2) 

So, 

/(a;o) = (l + x2')a;2'' 
= (l + a;2)^i 
= (1 + 3^)31 

= (10)3 = 30 = 13 mod 17 

The second, trickier equation comes from the same 
subsection. The algorithm to use is in Figure [10] on 
page m This is modified from Figure [7] on page |6l 
Here, 

9{^) = = 1 (4.3) 

There is nothing to multiply, so by convention this 
product will be set equal to one. As for x\. 

= {l + x^)x^ 
= (1 + 3^)3^ 

(4)9 = 36 = 2 mod 17 

Now, proceeding through the algorithm starting at 
h=0, that results in the case h = al. So the result 



becomes g{xo)^ or 1. Now, h is incremented, so h=l. 
This results in the case h — Ui. So g{xi) is added 
to the result, giving 1 + 2 ~ 3. This finishes the 
algorithm with the result f{xo V xi) = 3 mod 17. 

At some time, the equations should be set up so 
that Boolean algebra can eventually lead to an an- 
swer. Now is a good time, so the various results for 
ao — 2 and cq are shown in Figure [25] on the following 
page. The goal here is to eliminate the diagonal by 
taking combinations of equations, as explored in the 
previous section. Observe what happens if another 
set of equations is created, with oi = 3 and Ci = 1: 

ci • 1 = 1 • 1 =1 
ci • fli = 1 • 3 =3 
ci- 01^ = 1 -32 =9 

Proceeding as in Figure [23l on page[T7l the first dif- 
ferences of this equation are 3 — 1 = 2 and 9 — 3 = 6. 
The second difference is 6 — 2 = 4. Recalling that 
the second differences should add up to the modulus 
(17), a second difference of 17 — 4 = 13 is required. 
From Figure [25l it can be seen that setting co = 13 
in the first set of equations gives a second difference 
of 13. Thus, these two equations can now be com- 
bined succesfully. To see this, simply use the original 
values of a and c in each equation and add them to- 
gether. For example using ao — 2 and cq = 13 for the 
first equation, and ai = 3 and cq = 1 for the second 
equation, this gives results of combining equations: 

Co -I- ci = 13 + 1 =14 
co(ao) -}-ci(ai) = 9 + 3 =12 
co(ao2)+ci(ai2) = 1 + 9 =10 

Now corresponding addition equations can be con- 
structed by observing: 

14 ^=^ 14 + 0(-2) 
12 14 + l(-2) 
10 14 + 2(-2) 

A similar examination can be done to produce two 
more addition equations. Proceeding by increasing 
the previous value of a for each successive multiplica- 
tion equation, recall that the last multiplication equa- 
tion used fli = 3. So pick 02 = 4. Use C2 = 1 to start. 
This gives: 
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Figure 25: Multiplication Results with ao = 2 Modulo 17 
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C2 • 1 = 1 • 1 =1 
C2 • 02 = 1 • 4 =4 

C2 • = 1 • 4^ = 16 

This gives first differences of 4 — 1 = 3 and 16 — 4 = 
12, and a second difference of 12 — 3 = 9. So look for 
another second difference that sums 9 to a total of 
17. 17 — 9 = 8. Consulting Figure [55] again, it can be 
seen that using cq = 8 gives a second difference of 8. 
To check this and set up the values for the addition 
equations, use ag = 2, cq = 8, a2 = 4, and C2 = 1: 



This setup gives the information needed to set up 
a system of three equations in three unknowns. The 
next step involves preparing the clause polynomials 
for addition and multiplication. Once they are pre- 
pared, addition and multiplication can be used to set 
up the linear algebra system. 

First off, the function of ones should be calculated. 
This is easy, recall Equation 13. II from page fTTl 



Co + C2 = 8 + 1 =9 
co(ao) + C2(a2) = 16 + 4 =3 
co(ao^) + C2(a22) = 15 + 16 = 14 = -3 

Once again, observe that there is a difference of 
successive results by -6. For example 9 + (-6) = 3, /(I) = H (1 + 2; ) (4.4) 

and 3 + (—6) = —3. So the second set of addition '^"^ 
equations can be set up: 



9 9 + 0(-6) 
3 <^ 9 + l(-6) 
-3 <^ 9 + 2(-6) 

A third set of addition equations can now be set 
up, which will conclude the setup. Continuing suc- 
cessively, pick as — 5. With C3 = 1, this gives: 



Here the function of ones is: 



C3 ■ 1 = 1 • 1 
C3 ■ 03 = 1 • 5 



C3 • 02 



1-5" 



It gives first differences of 5 — 1 = 4 and 8 — 5 = 3, 
and a second difference of 3 — 4 = 16. 16 + 1 = 17, 
so a second difference of one is required for another 
set of equations. Once again using Figure [^H co = 1 
gives this equation for ao = 2. Thus: 

Co + C3 = 8 + 1 =2 
co(ao) + C3(a3) = 2 + 5 =7 
co(ao') + C3(a3') = 4 + 8 ee 12 

For this last set of equations we observe: 

2 +=^2 + 0(5) 
7+=^ 2 + 1(5) 
12 <^ 2 + 2(5) 



v-i 



/(i)=n (!+-') 

fe=0 

= {l+x'){l+x^) 
= (l + 3)(l + 9) 
= (4)(10)s6 



(4.5) 

(4.6) 
(4.7) 
(4.8) 
(4.9) 



Now that the function of ones is calculated, the 
multiplication equations can be calculated. 



to be completed later... 
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5 Finishing The Two Clause 
Example 

An important observation that is critical to the entire 
algorithm can be made at this point. Originally, the 
gist of the algorithm was to create 3 equations in 3 
unknowns. Unfortunately, the equations cannot be 
independent of one another, since we can only vary 2 
parameters of the 3 equations; Only the value for d, 
which corresponds with no clauses satisfied, and the 
corresponding increment can vary. Since only these 
two parameters can vary, we must somehow satisfy 
the 3 unknowns with only two equations. 

This is where the observation occurs. The three 
unknowns aren't completely independent. One of 

the 3 unknowns depends upon the other 2. 

In the case of addition, the three unknowns must 
add up to the function of ones. This is because all 
three unknowns, taken together, completely occupy 
all of the X coefficients, which is exactly what the 
function of ones represents. Similarly, in the case of 
multiplication, the off-diagonals must add up to the 
function of ones sqaured, minus the diagonal, which 
is the function of ones. 

Let's refine our model. First, we only need to come 
up with two equations. Back around pageHH and in 
Figure [Ml we came up with 3 variables; 6oj ^i, and 
&2- represents clauses satisfied, bi represents one 
clause satisfied, and 62 represents 2 clauses satisfied. 
We just realized that these three variables contain 
one dependent variable among them. Let's let 62 be 
the dependent variable, which makes the equation 



46o + 26i + O62 = 46o + 26i + 0(/(l) - 60 - ^'i) (5.1) 



Thus, we can write our old equation in 3 depen- 
dent variables as a new equation in 2 independent 
variables. So now we only need 2 equations, and 
we have them. We can now solve the system and 
determine the values associated with no clauses sat- 
isfied, one clause satisfied, and two clauses satisfied 
(/(I) — 60 — ^i)- This is just linear algebra, but care 
must be taken to ensure that only multiplications are 
performed instead of division, since we are working 
inside a finite field. 



6 Algorithm Finish (Basics) 



In all cases, the algorithm finishes the first portion 
with a value for n clauses satisfied. If this value is 
nonzero, we can conclude that the current Boolean 
equation can be satisfied. Otherwise, there is only a 
\-m-p2 chance that the current Boolean equation can 
be satisfied (where p2 is a probability picked ahead of 
time that will be discussed in greater detail shortly). 
But this doesn't return a certificate; that is, an as- 
signment of variables that satisfies the equation. If 
we think that the equation can be satisfied, we should 
return a certificate. 

To return a certificate, multiple occurences of the 
basic algorithm are run. This proceeds as follows. 
First, we assume that the equation can be satisfied; 
otherwise, simply return that it can't be satisfied and 
we are done. So the next step is to take the first 
value in the original equation and pick a value for it. 
We'll pick true, although we could pick false. Now we 
rewrite the original Boolean equation with this new 
value set as true (which is fairly common knowledge, 
and may be explained in a future version of this pa- 
per). Then we determine if the new equation can be 
satisfied. If it can, we proceed to repeat this process 
with more variables until a certificate is produced. If 
the equation can't be satisfied with the variable set 
as true, we try a new equation with the variable set 
as false. If this works, we again proceed on with more 
variables. Now there is a slight chance that neither 
equation seems to be satisfied. If this is the case, we 
can try again with another prime. What's important 
here is that there is no need to backtrack. 

This is a probabilistic algorithm with bounded er- 
ror. At the start of the algorithm, we must know the 
total probability of error, which we will call p^. Then 
we can calculate the maximum probability of error 
for each iteration, which is p2- Then we can deter- 
mine how many and what types of primes to use to 
give us p2 error at each step. All of this will be dis- 
cussed after the general n clause algorithm. For now 
we remark that we can continue to examine a par- 
ticular variable assignment only for so long, and then 
we can simply conclude that we've exceeded the error 
probability and conclude that we couldn't deduce a 
certificate; only satisfiability. 

The major point to take away here is that an n 
clause Boolean equation with V variables will even- 
tually require the algorithm for more clauses. Specif- 
ically, this is the n -\- 2V clause algorithm. 
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7 The n Clause Algorithm 

Knowing that the n clause Boolean equation will 

eventually require the n + 2V clause algorithm, it's 
best to anticipate this ahead of time. 

The actual general case algorithm works by using 
many versions of the 2 clause algorithm. Wc can be- 
gin to see how this can occur by observing that the 
2 clause algorithm returns a value for all clauses sat- 
isfied, which can be used as one of the four initial 
values in a second 2 clause algorithm. The second 2 
clause algorithm has the same general requirements 
as any 2 clause algorithm. It requires a prime, which 
we have. It requires the information for 2 clauses; 
two numbers for addition, and two numbers for mul- 
tiplication. Wc can get this by using 4 versions of the 
2 clause algorithm, which prepare another iteration 
of the 2 clause algorithm. Then we proceed again. In 
this fashion, a tree of 2 clause algorithms can be built 
up to solve an n clause system. 

7.1 Semi-optimized Version 

The version presented in this paper has only minor 
optimizations to enhance it; however, there will be 
much room left for further improvement. 

So far, the best way to optimize seems to be to 
focus on the main portion of the algorithm, which is 
the 2 clause algorithm. If wc briefly analyze the oc- 
curences of this, wc should know that it combines two 
clauses into one output. Each clause needs a portion 
defined for addition, and another for multiplication. 
So for every 2 clauses in, there are 4 total inputs. 
This then leads to a single output. So for n clauses, 
at every level the number of clauses is reduced by 
half. Similarly, the number of inputs is reduced to a 
fourth. We can designate the number of such levels 
as /. Now for n = 2' clauses there are I levels to the 
tree. Similarly, there are 4' inputs needed for all of 
the leaves. We note that the 4' inputs is equal to the 
2' clauses squared: 4' = (2')(2'), just as the inputs 
are the number of clauses squared. So we can con- 
clude that there are roughly total instances of the 
2 clause algorithm that occur for the n clause algo- 
rithm. We also know that the n clause algorithm will 
be run 0{V) times in the case of a satisfied Boolean 
equation in order to determine a certificate (a set of 
satisfying variables). So we can also conclude that 
the algorithm will run the 2 clause algorithm no more 
than 0{V){n + V)"^ times. 

Knowing that the current center of attention for 
the n clause algorithm is the repetition of the 2 clause 
algorithm, we can optimize the performance of the 2 
clause algorithm by making some precomputations. 



This is because wc can use; the; same values repeatedly 
for the 2 clause algorithm, with the only exceptions 
(outside of intermediate computations) being the 4 
inputs and single output. 

Knowing that the main algorithm will call the two 
clause algoithm no more than 0{V){n+V)'^ times, we 
can set all of the primes we use to be; approximately 
this value. Thus we set our primes that we use to be 
Q{V{n+V)^). All calculations can now be performed 
knowing the value of the prime ahead of time, and 
this will also help to ensure that the mulitplication 
inputs are rarely ever zero. In the case that they 
are, the algorithm will have to perform some addition 
steps to try to ensure that they become nonzero. 

Returning to the precalculations, we can set up the 
equations to proceed as quickly as possbilc. Calculate 
the appropriate a^'s and c;j's so that all two clause 
algorithms run with the same equations. In fact, all 
of the constants can be precalculated individually in 
0{V){n + V)"^ time by simply cycling through every 
possible natural (plus zero) up to the prime p and 
picking the appropriate value. 

So the precalculations all take 0{{V){n+V)'^) time. 
Note that the two clause algorithms, at this point, 
should take constant time. 

The algorithm is almost complete. Only one im- 
portant piece remains; the problem that arises when 
a zero is given as an input for multiplication. There 
is one good possibility that we can make use of. We 
use extra variables in our calculations. To do this, 
we simply perform clause calculations with double 
the number of variables in the equation (A different 
value could be used, but this seems to be a fairly use- 
ful amount). Now. when a zero comes up as an input 
for multiplication, we can simply add in a new equa- 
tion with the new variables. If the original equation 
is satisfiable, this should help to change the multipli- 
cation input. Otherwise, we will conclude that the 
equation is questionably unsatisfiable (It is unsatisfi- 
able to within the error probability 1/p). 

8 Runtime Analysis / 
Correctness 

As mentioned and explained in the previous section, 
most components take 0(V{n + V)'^) time. The error 
correction in the case of a zero input correction is 
seperate from the main algorithm, and obviously can 
be done in 0(V(n -|- V)"^) time. So for a single prime, 
the runtime is 0{V{n + V)'^). 

This is a nonrandomised, deterministic algorithm 
with bounded error. Each prime used gives an in- 
dividual error bound of l/Q{V{n + V)'^). Together, 
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P primes give approximately an error bound not cx- 9 Conclusion 
ceeding probability l/Q{V{n + V)"^)^. So the algo- 
rithm runs in time 0{P ■ V{n + V)"^) with mistaking 

satisfiable Boolean expressions as imsatisfiable with 
an approximate probablity l/ld{V{n + Vy)^ . 



I hope that this project presents sufficient evidence 
that P=NP. I've written the last few pages rather 
hastily, but hope to improve things soon. I'm start- 
ing on writing the code for this project, so that ev- 
erything can be put under better scrutiny. 
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